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Final  Technical  Report 

Objectives 

This  project  was  funded  through  an  AFOSR  AASERT  grant  and  supported  four  graduate 
students  to  work  with  Prof.  Harley  H.  Cudney.  The  overall  objective  of  the  grant  was  to 
develop  new  technologies  in  optimal  actuation  and  sensing  methods.  The  specific 
objective  of  the  grant  was  to  develop  a  working  relationship  with  an  Air  Force 
Laboratory  and  actively  perform  research  that  specifically  supported  their  mission. 

The  project  that  was  originally  identified  which  most  closely  supported  the  mission  of  the 
Structures  and  Dynamics  Branch  of  the  Air  Force  Phillips  Laboratory  at  Kirkland  Air 
Force  Base  was  to  develop  control  algorithms  to  test  on  the  space-flight  experiment 
ASTREX.  However,  another  requirement  that  the  project  must  satisfy  is  to  have  a 
research  topic  suitable  for  the  student  to  write  a  dissertation.  After  discussions  with  Dr. 
Alok  Das,  and  Capt.  Jeanine  Sullivan,  Ph.D.,  we  determined  the  ASTREX  project  did  not 
have  this  content.  In  subsequent  discussions,  we  identified  the  acoustic  control  of 
fairings  as  being  a  definite  technology  need  for  both  the  Air  Force  and  industry. 

To  fulfill  the  overall  objective  of  developing  new  technologies  for  optimal  actuation  and 
sensing,  we  investigated  developing  design  procedures  for  model  reference  adaptive 
control,  which  has  been  widely  adopted  in  the  structures  community  as  a  robust  control 
method.  However,  no  specific  design  procedures  exist.  Our  objective  was  to  develop  a 
design  procedure  and  assess  the  robustness  of  four  approaches  to  control:  a  classical  pole 
placement  method,  and  three  adaptive  control  approaches. 

Also  to  support  of  our  overall  objective,  an  adaptive  vibration  absorber  capable  of 
detecting  the  firequency  of  the  force  and  tuning  itself  automatically  to  that  particular 
frequency,  is  designed,  built,  and  tested. 

Status 

Design  methods  for  Model  Reference  Adaptive  Control:  Concluded 

The  Model  Reference  Adaptive  Control  Project  is  complete  and  summarized  below.  Mr. 
Jon  Hill  graduated  with  a  Master’s  Degree  in  Mechanical  Engineering  in  1995. 

Adaptive  Vibration  Absorbers:  Concluded 

An  adaptive  vibration  absorber  was  designed,  built,  and  tested.  A  model  was  developed 
and  this  model  was  used  to  develop  the  tuning  algorithm.  The  actively  controlled 
absorber  was  tested  and  shown  significant  performance  advantages  over  passive 
absorbers. 

Active  Acoustic  Control  of  Rocket  Payload  Fairings:  Continuing 


A  formulation  was  created  for  the  active  acoustic  control  problem.  We  consider 
cylinders  to  be  an  effective  model  of  a  rocket  fairing.  We  have  developed  an  analytical 
method  to  predict  the  noise  level  inside  the  cylinder  due  to  acoustic  excitation.  We  have 
also  created  a  cylinder  which  has  simply  supported  boundary  conditions.  Both  of  these 
steps  are  necessary  to  determine  the  control  authority  of  piezoelectric  actuators  for 
providing  effective  attenuation  of  the  acoustic  loads  due  to  the  rocket  being  launched. 

Accomplishments,  New  Findings,  and  Research  Highlights 

Design  methods  for  Model  Reference  Adaptive  Control 

In  this  study,  we  assessed  the  robustness  of  four  distinct  control  approaches:  pole 
placement;  the  command  generator  tracker  (CGT)  approach  to  model  reference  control; 
model  reference  adaptive  control  (MRAC);  and  MRAC  using  a  fixed  feedback  gain.  We 
use  a  second  order,  single-input  single-output  (SISO)  plant  to  examine  the  performance 
and  stability  of  each  method.  This  evaluation  spans  a  broad  range  of  design  goals  and 
uncertainty  in  models  of  the  plant. 

Pole  placement  and  CGT  designs  are  linear  and  relatively  easy  to  implement,  but  require 
explicit  knowledge  of  the  plant.  Although  MRAC  schemes  require  little  knowledge  of  the 
plant's  dynamic  characteristics,  such  algorithms  are  non-linear  and  involve  design 
variables  whose  effects  are  not  readily  apparent.  Currently,  there  are  no  general  design 
procedures  for  MRAC.  In  this  study,  we  propose  a  method  for  designing  an  MRAC 
controller  applied  to  a  second  order  SISO  plant.  This  method  does  not  require  the 
controller  to  be  tuned  for  different  closed-loop  performance  goals.  This  procedure  also 
creates  a  consistent  basis  for  comparing  the  robustness  of  all  four  algorithms. 

Pole  placement  and  the  CGT  control  perform  as  designed  if  the  plant  is  modeled 
correctly.  Under  this  circumstance,  the  adaptive  controllers  also  perform  at  levels 
equivalent  to  the  linear  algorithms.  However,  conditions  with  plant  modeling  error 
highlight  enormous  differences  among  the  four  algorithms.  Pole  placement  suffers  the 
largest  response  error  and  for  extreme  testing  conditions,  instability.  The  CGT  controller 
exhibits  better  performance  than  pole  placement  and  remains  stable  over  all  testing 
variables.  MRAC  maintains  a  high  performance  level  under  severe  testing  conditions. 
MRAC  requires  minimal  plant  knowledge  to  guarantee  stability  and  good  performance. 

Our  model  reference  adaptive  control  design  procedure  yielded  excellent  results  in 
simulation.  Considerations  of  representative  model  trajectories,  representative  output 
error,  sampling  time,  and  transient  response  time  factored  into  our  weighting  matrix 
selection.  We  used  Lyapunov  stability  to  results  to  determine  a  feedthrough  term, 
minimizing  the  effects  of  plant  modeling  error.  In  the  presence  of  plant  modeling  error, 
these  design  considerations  were  sufficient  to  produce  responses  well  correlated  with  the 
reference  model’s.  The  model  reference  adaptive  controller  performed  particularly  well 
over  long  transient  periods,  where  the  adaptive  gains  could  accurately  converge  to  values 
minimizing  output  error. 


A  major  benefit  of  our  design  procedure  is  the  absence  of  user-controller  interaction.  This 
procedure  only  requires  an  approximate  model  of  the  plant  and  a  specified  reference 
model.  Although  our  study  assumes  a  very  simple  system,  where  the  reference  model 
and  plant  are  both  second  order,  this  design  method  is  innovative  and  may  be  applied  to  a 
larger  class  of  systems. 

In  the  comparative  analysis  of  robustness,  we  found  that  generally,  adaptive  control 
exceeds  the  ability  of  linear  controllers  to  track  the  response  of  the  reference  model. 
Although  the  linear  model  reference  control  method  yielded  lower  performance  indexes 
periodically,  the  adaptive  methods  begin  to  track  the  model  output  in  fewer  time  steps 
than  the  linear  methods.  This  attribute  of  adaptive  control  is  advantageous,  but  not 
sufficiently  recognized  by  our  quantitative  evaluation. 

We  raised  two  issues  worthy  of  further  investigation.  First,  we  should  examine  the 
effects  of  signal  noise  for  a  particular  choice  of  representative  error.  Our  selection  was 
very  low,  resulting  in  a  relatively  high  feedback  weighting  value.  This  high 
computational  value  did  not  pose  a  problem  in  our  “sterile”  testing  environment. 
However,  it  is  possible  to  achieve  equitable  performance  for  lower  feedback  weighting 
values. 

The  two  types  of  MRAC  schemes  produced  similar  results.  We  hypothesized  that 
MRAC,  using  a  fixed  feedback  gain,  would  have  a  lower  sensitivity  to  noise  than  the 
fully  adaptive  version.  The  fully  adaptive  MRAC  performed  marginally  better  than  its 
counterpart  in  our  study.  We  should  reexamine  these  methods  under  more  realistic 
testing  conditions. 

We  did  not  attempt  control  for  an  input  to  the  reference  model.  Our  design  procedure  can 
be  easily  modified  to  account  for  this.  Using  a  reference  model  command,  the  plant  has  a 
non-zero  steady  state.  This  gives  the  designer  a  perspective  to  better  evaluate  and  design 
the  ratio  between  proportional  and  integral  weighting  matrices.  A  non-zero  steady  state 
also  accentuates  the  benefits  of  adaptive  control  over  linear  control. 

Adaptive  Vibration  Absorbers 

An  adaptive  vibration  absorber  capable  of  detecting  the  fi’equency  of  the  force  and  tuning 
itself  automatically  to  that  particular  frequency,  is  designed,  modeled,  and  tested.  The 
design  stage  begins  with  a  criteria  list  that  evaluates  various  design  concepts.  After  a 
concept  is  chosen,  the  concept  is  developed  using  vibration  theory  and  beam  equations. 
The  principal  system  is  designed  first  to  obtain  a  certain  natural  frequency.  The  absorber 
is  designed  next  so  that  its  range  of  frequencies  includes  the  natural  frequency  of  the 
principal  system.  The  design  consists  of  a  cantilever  beam  as  the  principal  system,  and 
an  assembly  of  three  rods  and  a  motor  as  the  vibration  absorber. 

Preliminary  component  testing  was  performed  to  develop  a  model  of  the  principal 
system,  absorber  system,  and  the  combined  (principal  and  absorber)  system.  The  models 
predict  the  natural  frequencies  of  the  systems,  and  in  the  case  for  the  combined  system, 
the  model  predicts  the  anti-resonant  and  resonant  frequencies.  Experiments  of  the 


combined  system  were  performed  to  evaluate  the  effectiveness  of  the  absorber. 
Furthermore,  experimental  results  for  all  three  systems  were  compare  to  the  results  from 
the  theoretical  model. 

The  combined  system  evaluation  shows  that  the  absorber  works  effectively.  It  is  able  to 
include  the  natural  frequency  of  the  principal  system  in  its  range  of  frequencies.  The 
model  predicts  the  anti-resonant  frequency  within  a  mean  of  0.3  %  error.  Furthermore, 
the  model  predicts  the  two  resonant  frequencies  (coi  and  02)  within  a  mean  of  0.6  %  error 
and  0.6  %  error,  respectively. 

A  control  algorithm  is  developed  and  applied  to  this  tunable  vibration  absorber.  The 
adaptive  vibration  absorber  is  capable  of  detecting  the  frequency  of  the  driving  force  and 
tuning  itself  automatically  to  that  frequency.  The  primary  structure  was  previously 
designed  to  obtain  a  certain  natural  frequency,  and  the  absorber  was  designed  such  that 
its  range  of  frequencies  included  the  natural  frequency  of  the  primary  structure.  The 
primary  structure  consists  of  a  cantilevered  beam  with  the  absorber  attachment  hardware, 
and  the  vibration  absorber  assembly  consists  of  three  rods  and  a  stepper  motor. 

The  control  algorithm  uses  a  look-up  table  and  a  gradient  search  to  optimize  the 
effectiveness  of  the  absorber  to  reduce  the  vibrations  on  the  primary  structure.  The  look¬ 
up  table  uses  an  equation  which  is  based  on  experimental  data  to  transform  a  given 
voltage  input  (directly  proportional  to  the  forcing  frequency)  into  an  output  command 
necessary  to  adjust  the  natural  frequency  of  the  absorber.  One  the  input  voltage  reaches 
steady  state,  the  gradient  search  routine  adjusts  the  natural  frequency  of  the  absorber  to 
ensure  the  absorber  is  tuned  to  the  optimal  frequency  that  minimizes  the  vibration 
amplitude  of  the  primary  structure. 

The  primary  structure  with  the  adaptive  absorber  offers  significant  reductions  in  the 
vibration  amplitudes  of  the  primary  structure,  as  compared  to  the  primary  structure  with 
no  absorber,  and  to  the  primary  structure  with  a  passive  absorber. 

Active  Acoustic  Control  of  Rocket  Payload  Fairings 

Future  launch  vehicle  payload  fairings  will  be  manufactured  from  advanced  lightweight 
composite  materials.  The  loss  of  distributed  mass  causes  a  significant  increase  in  the 
internal  acoustic  environment,  causing  a  severe  threat  to  the  payload.  During  launch  the 
payload  fairing  is  excited  by  a  severe  acoustic  loading  provided  from  the  bum  of  the 
rocket  engines.  The  acoustic  noise  is  transmitted  through  the  fairing  and  creates  an 
internal  SPL  ranging  from  120  to  140  dB  that  can  damage  the  payload.  There  remains  a 
great  deal  of  work  in  understanding  the  Physics  behind  the  acoustical  problem  and  the 
interaction  between  the  external  and  internal  SPL’s  and  the  payload  fairing  skin.  Current 
state  of  the  art  active  control  technology  has  been  considered  to  tackle  the  fairing 
acoustical  problem.  The  actuator  types  evaluated  include  acoustical  noise  sources,  point 
forces,  magnetostrictive  actuators,  smart  materials,  electrostrictive  and  piezoelectric 
actuators.  Currently  the  best  candidate  to  solve  the  problem  is  the  piezoelectric  actuator. 
However,  due  to  the  broadband  and  intense  nature  of  the  acoustic  environment,  it  remains 
uncertain  if  piezoelectric  actuators  will  have  the  authority  to  solve  the  problem. 


To  help  determine  the  acoustic  control  authority  of  piezoelectric  actuators  mounted  on  a 
rocket  fairing,  the  internal  acoustic  response  created  by  the  actuators  needs  to  be 
determined.  The  internal  acoustic  response  of  a  closed  simply-supported  cylinder 
actuated  by  piezoelectric  actuators  can  be  conveniently  determined  by  using  the  Kirchoff- 
Helmholtz  integral  and  the  cylinder’s  spatial  surface  velocity.  The  surface  velocity  is 
determined  using  an  impedance  model  of  an  actuator  bonded  to  the  surface  of  a  cylinder. 
The  structural  analysis  is  compared  to  an  experiment  and  the  results  are  in  agreement. 
The  acoustic  formulation  is  described  along  with  a  description  of  future  work.  The  work 
presented  provides  the  initial  steps  required  in  determining  the  acoustic  authority  and 
feasibility  of  using  piezoelectric  actuators  on  rocket  fairings. 

The  structural  model  we  developed  needs  to  be  further  verified  for  a  broader  firequency 
range  with  more  experimental  information  about  the  absolute  spatial  displacement  of  the 
cylinder.  This  can  be  accomplished  by  using  scanning  laser  vibrometry.  The  model  also 
needs  to  be  verified  for  out  of  phase  actuation  and  more  attention  needs  to  be  given  to  the 
actual  damping  values  on  a  mode  for  mode  basis.  The  results  from  the  impedance  model 
will  ultimately  be  used  to  calculate  the  acoustic  response  within  the  cylinder.  Once  the 
acoustic  model  is  verified,  the  combined  structural/analytical  model  will  be  used  to 
determine  the  authority  of  PZT  actuators  on  a  cylinder  with  properties  similar  to  a  rocket 
payload  fairing. 

An  experiment  of  a  cylinder  with  simply-supported  boundary  conditions  has  been 
conducted  and  was  shown  to  exhibit  d5mamic  behavior  very  similar  to  what  the  theory 
predicts.  An  analytical  impedance  based  model  of  a  PZT  actuating  on  a  SS  cylinder  was 
also  described  and  compared  with  the  experimental  results.  The  results  indicate  that  the 
impedance  based  model  can  predict  the  overall  resonant  and  modal  behavior  but,  more 
importantly  for  acoustical  applications,  is  able  to  determine  the  magnitude  of  the 
vibration.  An  acoustical  formulation  of  the  pressure  field  within  the  cylinder  has  also 
been  described  based  on  the  Kirchoff-Helmholtz  integral.  The  work  presented  provides 
the  initial  steps  required  in  determining  the  acoustic  authority  and  feasibility  of  using 
piezoelectric  actuators  on  rocket  fairings. 

In  many  aerospace  and  strucUiral  applications,  idealized  cylinders  are  used  to  approximate  more 
complex  structures  such  as  airplane  fuselages  and  rocket  payload  fairings.  Many  authors  have 
created  sophisticated  models  of  the  structural  response  of  cylinders  excited  by  a  variety  of 
actuators.  However  seldom  are  these  models  verified  experimentally.  None  of  the  models 
created  to  describe  the  response  of  a  simply-supported  (SS)  cylinder  have  ever  been  verified.  In 
this  work  several  different  boundary  conditions  are  created  to  approximate  an  ideal  SS  cylinder. 
The  different  boundary  conditions  tested  are  described  and  compared  with  finite  element  analysis 
(FEA)  and  with  an  impedance  based  analytical  model  for  a  piezoelectric  actuator  (PZT)  exciting  a 
SS  cylinder.  The  work  presented  describes  the  steps  taken  in  creating  a  physical  boundary 
condition  which  approximates  the  ideal  boundary  condition  for  a  SS  cylinder.  The  results 
indicate  that  the  created  physical  boundary  condition  resembles  the  ideal  SS  boundary  condition. 
The  structural  response  computed  by  the  impedance  model  is  also  shown  to  agree  with  the 
experiment  and  with  the  FEA.  The  results  can  be  used  to  create  other  experimental  SS  cylinders 
and  will  be  used  at  Virginia  Tech  to  determine  the  feasibility  of  applying  piezoelectric  actuators 
to  control  rocket  payload  fairings  vibrations. 


The  work  presented  describes  the  steps  taken  in  creating  a  physical  boundary  condition  which 
approximates  the  ideal  boundary  condition  for  a  SS  cylinder.  The  results  indicate  that  the  created 
physical  boundary  condition  resembles  the  ideal  SS  boundary  condition  but  is  not  perfect.  It  does 
however  have  very  similar  vibrating  properties  and  modal  distribution.  The  analytical  model 
predicts  the  expected  physical  changes  to  the  response  as  the  actuator  location  is  shifted.  The 
structural  response  computed  by  the  impedance  model  is  also  shown  to  be  similar  to  the 
experimental  response  and  matches  the  natural  frequencies  determined  using  FEA.  The  results 
can  be  used  to  create  other  experimental  SS  cylinders  and  will  be  used  at  Virginia  Tech  to 
determine  the  feasibility  of  applying  piezoelectric  actuators  to  control  rocket  payload  fairings 
vibrations. 
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Abstract 

Future  launch  vehicle  payload  Imrings  will  be 
manufactured  from  advanced  lightweight  composite 
materials.  The  loss  of  distributed  mass  causes  a 
significant  increase  in  the  internal  acoustic 
environment,  causing  a  severe  threat  to  the  payload. 
Using  piezoelectric  actuators  to  control  the  fairing 
vibration  and  the  internal  acoustic  environment  has 
recently  been  proposed.  The  control  authority  of  these 
actuators  for  this  problem  has  not  yet  been  determined. 
To  help  determine  the  acoustic  control  authority  of 
piezoelectric  actuators  mounted  on  a  rocket  fairing  the 
internal  acoustic  response  created  by  the  actuators 
needs  to  be  determined.  The  internal  acoustic  response 
of  a  closed  simply-supported  cylinder  actuated  by 
piezoelectric  actuators  can  be  conveniently  determined 
by  using  the  Kirchoff-Helmholtz  inte^  and  the 
cylinder’s  spatial  surface  velocity.  The  surface  velocity 
is  determine  using  an  impedance  model  of  an  actuator 
bonded  to  the  surface  of  a  cylinder.  The  stmctural 
analysis  is  compared  to  an  experiment  and  the  results 
are  in  agreement  The  acoustic  formulation  is  described 
along  with  a  description  of  future  work.  The  work 
presented  provides  the  initial  steps  required  in 
determining  the  acoustic  authority  and  feasibility  of 
using  piezoelectric  actuators  on  rocket  fairings. 

Introduction 

Traditionally,  payload  fairings  have  been  constructed 
from  lightweight  aluminum  alloys.  However,  designers 
are  currently  developing  composite  fairings  which  are 
much  easier  to  manufacture  while  also  being  lighter  and 
stiflfer  than  existing  aluminum  fairings. 

Reducing  the  acoustic  noise  that  the  payload  (i.e. 
satellite)  receives  during  launch  has  always  been  an 
important  design  constraint  An  excess  of  acoustic 
noise  transmitted  through  the  fairing  poses  a  potential 
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threat  to  the  life  of  the  payload.  Since  future  fairings 
will  be  constructed  of  lighter  materials  having  a  lower 
transmission  loss  (TL),  the  acoustic  threat  is  even  more 
severe.  The  most  recent  research  for  the  problem  has 
implemented  passive  techniques  in  the  form  of 
improved  blankets  and  ^mamic  vibration  absorbers 
mounted  on  the  payload’"* **. 

Active  control  technology  is  currently  being 
investigated  as  a  potential  solution  for  the  fairing 
acoustic  problem  in  the  low  frequency  range  (50400 
Hz)  where  blankets  are  less  effective.  A  review  of 
ciurent  active  control  actuators  has  been  recently 
performed  along  with  a  description  of  the  acoustic 
fields  needing  control®.  The  use  of  piezoelectric  (PZT) 
actuators  to  actively  control  fairing  stmctural  vibrations 
and  the  internal  acoustic  environment  has  recently  been 
proposed.  However,  the  authority  of  these  actuators  for 
this  problem  has  not  yet  been  determined.  To  help 
determine  the  acoustic  control  authority  of  PZT 
actuators  mounted  on  fairing  like  stmctures,  the  internal 
acoustic  response  created  by  the  actuators  needs  to  be 
determined.  For  a  clos^  cylinder,  this  can  be 
accomplished  by  using  the  Kirchoff-Helmholtz  (K-H) 
integral  once  the  spatial  stmctural  response  is 
determined.  A  similar  analysis  has  been  performed,  but 
only  implemented  harmonic  point  load  excitation^ 

Presented  within  this  paper  is  a  description  of  the 
theoretical  analysis  us^  to  determine  die  stmctural 
response  of  a  simply-supported  (SS)  cylinder  via  an 
impedance  based  model  for  the  PZT  actuator.  The 
theoretical  solution  of  the  internal  acoustic  field  using 
the  stmctural  response  is  also  discussed  along  with 
some  advantages  of  using  this  technique.  An 
experiment  is  performed  to  validate  the  stmctural 
m^el  and  the  empirical  and  analytical  results  are 
compared.  To  the  author’s  knowledge  an 
analytical/experimental  comparison  of  a  SS  cylinder 
has  not  been  performed.  Finely  the  major  findings  and 
future  work  are  presented. 


Theoretical  Development 
PZT  Actuator  Model  and  Structural  Response 
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Theoretical  Development 

P7.T  Actuator  Model  and  Structural  Response 

Several  models  for  PZT  actuator  patches  attached  to 
cylinders  have  been  presented.  Initially  plate  models 
were  adapted  to  (flinders  via  the  Donnel-Mushtari  and 
Love  shell  equations*’’.  Later  models  were  based  on 
shell  governing  equations'®"**.  All  of  these  models  are 
based  on  force  inputs  calculated  at  zero  frequency 
excitation.  None  of  them  take  into  consideration  the 
changing  impedance  of  the  structure  as  a  function  of 
fiequency.  This  will  lead  to  an  incorrect  response  if 
the  structural  impedance  changes  with  frequency. 
Impedance  models  of  PZT  actuators  have  been 
previously  developed  for  beams,  plates,  rings  and 
shells'*”.  It  has  been  shown  that  neglecting  the 
frequency  varying  structural  impedance  can  lead  to  an 
incorrect  response.  The  structural  response  of  a  SS 
cylinder  is  determined  using  the  impedance  model 
developed  by  Lalande'®.  A  general  overview  of  this 
model  will  now  be  presented. 

The  coordinate  system  used  for  the  cylinder  and 
actuator  is  shown  in  Fig.  1.  The  model  assumes  two 
rectangular  PZT  actuators  are  attached  to  the  cylinder, 
one  on  the  outer  surface  and  another  collocated  on  the 
inner  surface.  The  resultant  forces  (in  phase  actuation) 
and  moments  (out  of  phase  actuation)  are  assumed  to 
act  at  the  cylinder’s  midplane  as  shown  in  Fig.  2. 


M,=(t,^t,)F,  (Id) 


Figure  2.  Transfer  of  Actuator  Forces  to 
Equivalent  Mid-Plane  Loading 

where  t  and  /  are  the  actuator  and  shell  thickness, 

<M  S 

respectively.  The  actuator  forces  and  Fq  need  to  be 

solved  by  using  the  PZT  actuator  constitutive 
equations,  actuator  impedance  and  the  structural 
impedance  or  admittance  (see  Fig.  3).  The  structural 
impedance  along  the  actuator  edge  is  related  to  the 
structural  admittance  by; 

Zxx 


Uxx 

H„ 


-1 


(2) 


Figure  1.  Cylinder  with  Mounted  Actuators 

The  relationship  between  the  actuator  output  force  and 
the  equivalent  midplane  force  or  moment  in  the  axial 
or  tangential  direction  is  given  by; 

inphase; 

=  2F,  (la) 

Ng  =  IFg  (lb) 

out  of  phase; 


where  Zand  represent  the  structural  impedance 
and  admittance,  respectively.  The  structural 
admittance  will  be  discussed  later.  The  input 
impedance  (shorted)  of  the  PZT  actuators  are  given  by. 


direct  terms; 

““  j<»tan(xfJ 
ia  tan(KR„ ) 


(3a) 

(3b) 


cross  terms; 


^axe  ~ 


ieo  tanitdj 


(3c) 


iG)tan(KRJ 


(3d) 


is  the  actuator  complex  modulus  of 


elasticity  with  no  applied  electric  field,  is  the 
actuator  angular  length,  is  the  actuator  axial  length. 
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<D  represents  the  angular  frequency,  /  is  an  imaginaiy 
number,  and  K  is  given  by: 

(4) 

where,  is  the  density  of  the  actuator. 


Figure  3.  Interaction  Between  the  PZT  and  the 
Cylinder  via  the  Structural  Impedance 

Combining  the  actuator  impedances,  Eq.  3,  with  the 
PZT  constitutive  equations  and  ^plying  boundary 
conditions,  a  relationship  between  the  dynamic 
actuator  force  and  the  stnictuial  impedance  (Eq.  2) 
results*’*’: 


F,  =-^(^sin«)Z„  +Csin(icRJZ,,)e'"  (5a) 

a 

F«  =  sm« )Za,  +  Csin(/cR, (5b) 

a 

where,  A  and  C  are  calculated  from. 


/£:cos(idj|^ 

(z . 

/ccos(Ad,)  — ^ 


1- 

Z 


ixO 

-  V. 


iccos(/cR,) 


-  V. 


a:cos(kR^) 


1-v. 


(6) 


where,  V ^ ,  c/j,  and  afjj  are  the  actuator  Poisson  ratio 
and  PZT  constants  and  Ej.  is  the  applied  electric  field. 

By  using  Eqs.  1  and  5,  the  actuator  output  forces  and 
moments  can  be  used  to  determine  the  structural 


response  assuming  the  structural  impedance  is  known. 
The  computation  of  the  structural  response  and 
structural  impedance  shall  be  presented  next. 


To  determine  the  structural  response  of  the  SS  cylinder 
caused  Ity  the  actuation  of  the  PZT  actuator,  a  modal 
expansion  analysis  is  performed.  The  equations  that 
describe  the  motion  of  the  cylinder  arc  based  on  Love’s 
equations  for  shell  structures.  The  general  equations 
have  been  adapted  to  incorporate  loading  caused  by  the 
actuation  of  line  moments  or  forces^®.  Subsequently, 
the  equations  are  developed  to  incorporate  the  loading 
generated  by  PZT  actuators* 


For  a  cylinder,  the  natural  orthogonal  modes  can  be 
used  to  determine  the  forced  response  using  an  infinite 
series: 

=  (7) 

where,  the  subscript  /  corresponds  to  the  associated 
displacement  (i.e.  /  =  3  represents  transverse 

displacement),  p^O)  represents  the  modal 
participation  factor  and  U^(x,ff)is  the  spatial  mode 
shape.  For  a  SS  cylinder  the  mode  shape  is  given  by’®: 


U,i(x,  0)  =  A^  cos(ax)cos(n0) 

(8a) 

Ujt(x,  9)  =  B^  sin(ca)sin(M^) 

(8b) 

C/j*  (x,  =  sin(aa:)  cos(n^) 

(8c) 

mn 

(8d) 

where  subscripts  m  and  n  represent  the  modal  indices 
and  subscript/?  refers  to  one  of  the  three  principal 
coordinates  (1,2,3).  After  performing  the  standard 
modal  expansion  procedure  to  the  differential 
equations  (multiplying  by  the  mode  shape,  integrating 
over  the  domain,  and  applying  orthogonality 
conditions),  the  modal  participation  &aor  can  be 
determined  by  solving  an  ordinary  differential  equation 
with  a  redefined  forcing  function: 

A(0  +  6)’/)*(0  =  F,c-  (9) 


where,  is  the  a  modified  forcing  function 
determined  by  the  actuator  forcing  terms  in  the 
differential  equation.  The  modified  forcing  function  is 
defined  in  reference  20  and  not  included  here  for 
brevity.  The  solution  to  the  modal  participation  factor 
can  be  given  in  the  form: 


(10) 
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Sg  =sin«^,  -sin«^j 

Cg  =COS«^,  -C0S«^J 

=  cos  CCCi  -  cos  CCCj 

where,  N’^  is  the  mode  normalization  constant: 
inphase: 


XT-  BIk 

out  of  phase: 


A. 


V  ^ imp 


+ 


B  ^  ^ 

^  I 

^imp  J 


«  ?!:  0,  /W  0 


XT-  BIk 


B^^ 

I 


« 5* 0,  m^O 


The  ratios  jC^  can  be  found  in 

reference  20.  Note  that  reference  18  and  19  have  some 
typographical  errors  that  have  been  corrected  here. 
Finally  the  structural  impedance  needs  to  be  described. 
For  two  dimensional  structures,  the  admittance  matrix 
is  given  by’®: 


/ 

• 

\ 

• 

aw 

2u 

irphase  GC 

V 

outphaseJ 

itphase 


ac 


=  -{h„f^+hm 


jr=xi 


(12a) 


=- 
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J  00  00 

ZEE 


/?(^2  —  )  p=i  „-i 


Xi. 


(  A  ^ 


V  ^  imp  J 


where  0^  refers  to  the  actuator  center.  The  structural 
admittance  (Eq.  13)  can  now  be  used  with  Eq.  2  to 
determine  the  structural  impedance,  and  therefore  the 
actuator  forces  can  be  determined  and  substituted  into 
Eq.  11  to  determine  the  transverse  displacement  of  the 
cylinder.  It  is  the  transverse  velocity  component  of  the 
cylinder  that  generates  the  acoustic  field,  and  is  what  is 
of  interest  in  this  work.  The  calculation  of  the  internal 
acoustic  response  shall  be  addressed  next. 


Acoustic  Analysis 

One  way  to  compute  the  acoustic  field  produced  by  the 
vibration  of  a  structure  is  to  use  the  Kirchoff-Helmholtz 
(K-H)  integral.  If  the  structural  surface  vibration  is  the 
only  source  of  acoustic  energy,  the  K-H  integral  is 
given  by^': 


(14) 


2v 


inphase  ^  RdO 


\ 

-(H 

outphaseJ 

0=02 

inphase 


dw 


Rm 


ouphase  J 


=  <H^gF,-^HggFg)  (  1  2b) 


e=0x 


where, 

^  2 

Based  on  the  structural  admittance  defined  in  Eq.  12, 
the  direct  structural  admittances  are  determined  by 
dividing  the  axial  or  angular  velocity  by  the  appropriate 
force,  or  Fg.  For  brevity,  only  77^  is  presented. 
The  remainder  of  the  admittances  can  be  found  in 
reference  18  or  19. 


where,  p,r,r^  ,n,p^  ,w^  ,&G,  represent  the 
pressure,  field  point  location,  surface  point  location, 
outward  normal  to  the  local  field  velocity,  density  of 
the  fluid  medium  (air),  local  normal  sur&ce  velocity, 
and  acoustic  Green’s  fimction,  respectively.  The 
Green’s  fimction  must  satisfy  the  acoustic  wave 
equation  and  the  far-field  radiation  conditions.  The  K- 
H  integral  can  then  be  directly  integrated  to  solve  for 
the  acoustic  field  generated  by  a  structure  as  long  as  the 
surface  velocity  and  the  surface  pressure  is  known. 
Generally  the  spatial  siuface  velocity  of  a  structure  can 
be  measured  or  determined  via  finite  element  analysis 
or  using  analytical  techniques.  However,  the  spatial 
surface  pressure  is  not  easy  to  measure  and  is  difficult 
to  calculate. 

Presently,  acoustic  boundary  element  analysis  divides 
the  surface  into  N  elements  and  uses  the  integral  to 
calculate  the  pressure  at  the  surface.  N  simultaneous 


4 

American  Institute  of  Aeronautics  and  Astronautics 


equations  need  to  be  solved  and  this  leads  to  some 
mathematical  diEBculties.  Using  the  firee-field  Green’s 
function,  the  integral  contains  singularities  when 
r  =  r,.  For  the  radiation  problem  the  formulation  also 
fails  to  provide  a  unique  solution  to  the  surface  pressure 
at  interior  acoustic  eigen  frequencies  of  the  stmctural 
surface.  Several  methods  have  been  presented  to  solve 
these  difficulties  for  the  radiation  problem^^  ^  The 
solution  of  the  internal  problem  however,  has  not 
received  as  much  attention^’^^'^. 

For  the  internal  acoustic  problem,  a  convenient  form  of 
the  Green’s  function  can  facilitate  in  solving  for  the 
internal  pressure  field  within  a  closed  space.  By 
choosing  the  Green’s  function  such  that, 

=  0  (15) 

The  first  integral  in  Eq.  14  is  eliminated  and  the 
calculation  of  the  surface  pressure  is  no  longer  required. 
Therefore  the  pressure  field  can  be  solved  by  simply 
knowing  the  normal  surface  velocity  of  the  structure. 
The  Green’s  function  that  satisfies  Eq.  15  is  expressed 
in  terms  of  the  acoustical  modes  of  the  cavity  and  is 
given  by^': 


G(F,A;,m)  =  J 

ijk 


¥  ijk(r)¥  ijkil) 


where. 


c(F)  = 


1  within  V 

2  on  the  boundary  of  V 

K.=\¥]ArW 


(16a) 


(16b) 


(16c) 


where  A  and  subscript  ijk  are  the  internal 


acoustic  mode  shape,  wavenumber,  modal 
normalization  constant  and  modal  indices,  respectively. 
For  complicated  geometries,  the  acoustic  mode  shape  is 
difficult  to  compute  and  so  direct  calculation  of  Eq.  14 
is  generally  not  possible.  However,  for  a  finite  closed 
cylinder  the  acoustic  mode  shapes  and  eigen  values  are 
known.  The  acoustic  mode  shape  is  given  by^’ : 


¥  ijk  =  Jj  (^jk  rlR)  cosiiTtxIl)  co%(j9)  (17) 

where,  J represents  a  j*  order  Bessel  function  and 

X  are  the  associated  roots  of  the  Bessel  function.  The 
modal  normalization  constant  is  the  mean  square  value 
of  the  acoustic  mode  shape  multiplied  by  the  volume  of 
the  cylinder.  For  a  cylindrical  duct  the  mean  square 
value  of  2-D  mode  shape  is  given  for  it’s  cross 
section^^.  This  expression  can  be  multiplied  by  the 


cylinder’s  axial  mode  shape  and  integrated  over  the 
length  of  the  cylinder  to  yield: 


where. 


e,ej 


e, ,  Sj 


1  /  =  0,y  =  0 

2  />o,y>o 


(18) 


and  is  related  to  the  roots  of  the  Bessel  function. 

All  the  terms  in  the  K-H  integral  (Eq.  14)  have  now 
been  defined  and  the  pressure  within  the  cylinder  can 
be  determined  directly.  The  only  difficulty  in  this 
formulation  is  the  singularity  of  the  Green’s  function 
when  the  wavenumber  coincides  with  one  of  the 
acoustic  eigen  fioquencies  of  the  cylinder.  This  can  be 
addressed  by  including  a  small  complex  component  in 
the  acoustic  wave  speed“. 


Advantaues/Disadvantages 

The  analysis  presented  in  this  worilc  has  several  distinct 
advantages  over  existing  techniques  as  well  as  some 
drawbacks.  The  structural  response  is  based  on  a  modal 
expansion  formulation  using  the  impedance  model  for 
PZJ  actuators.  An  impedance  based  model  of  the 
actuators  is  a  significant  improvement  over  existing 
models  which  do  not  consider  the  changing  structural 
impedance  of  the  cylinder  with  floquency.  The 
analytical  formulation  also  is  much  more  flexible  to 
changes  in  the  parameters  of  the  analysis  (i.e.  actuator 
size,  position  and  loading)  as  compared  to  performing 
an  equivalent  finite  element  analysis  (FEA).  As  with 
all  structural  analysis  based  on  modal  expansion 
(including  FEA),  the  choice  of  damping  affects  the 
absolute  response  of  the  cylinder.  l>amping  values  are 
very  difficult  predict  and  have  been  shown  to  vary  firom 
one  mode  to  ffie  next.  The  analysis  presented  assumes 
a  constant  damping  ratio  and  this  constitutes  an 
approximatioa  For  the  acoustic  response,  the  direct 
integration  of  the  K-H  integral  reduces  the  order  of  the 
problem  as  compared  with  using  FEA  The  analysis 
also  circumvents  calculating  the  surface  pressure  on  the 
inside  of  the  cylinder  to  determine  internal  acoustic 
response  as  is  required  by  boundary  element  methods. 
This  eliminates  the  singular  integral  caused  when  the 
field  point  coincides  with  the  surface  point  However, 
due  to  expanding  the  Green’s  function  in  terms  of 
acoustic  mode  shapes,  a  singularity  exists  when  the 
excitation  fiuquency  coincides  with  an  eigen  frequency 
of  the  enclosed  cavity.  The  singularity  is  a 
mathematical  deficiency  and  does  not  exist  physically, 
but  can  be  handled  by  introducing  an  imaginary 
component  to  the  acoustic  wave  speed. 
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Experiment/Anaivsis 


The  experimental  setup  used  to  model  a  SS  cylinder  is 
shown  in  Fig.  4. 


♦  Note;  All  values  should  be  converted  to  SI  units 

Table  1.  Cylinder  and  Actuator  Properties 


Loss  Factor 

0.005 

0.001 

16 

1.5 

■ssaiBMda 

10 

- 

0.25 

0.095 

i^hsssibs^hi 

- 

1.24 

- 

50  Vrms 

d32  (m/V) 

- 

-166xl0'* 

Figure  4.  Experimental  SS  Cylinder  with  Collocated 
PZT  Actuators  Located  at  3f=16/3”,  6  —  0° 

The  structure  consists  of  an  aluminum  cylinder  with 
two  PZT  actuators  (G1195)  collocated  at  x=16/3”, 

^  =  0° ,  driven  in  phase.  The  boundary  conditions  are 
created  by  two  thin  aluminum  annular  shims  (0.020” 
thick)  and  two  thick  aluminum  circular  disks  (.5”  thick, 
8.75”  diam.).  The  shims  are  screwed  to  the  cylinder 
ends  and  bolted  (w/epoxy)  to  the  disks.  Table  1 
describes  the  properties  of  the  instrumented  cylinder.  A 
roving  accelerometer  test  was  performed  on  the 
cylinder  while  the  PZT’s  were  driven  (in  phase)  at  50 
Vrms  from  750  to  2250  Hz.  A  Kistler  86161A500 
accelerometer  was  used  to  measure  the  cylinder 
response  at  x=l/4  and  at  x=l/2  for  various  angular 
positions  (0,  15,  30,...  165,  180  deg.)  which  were 
mirrored  about  180  deg.  A  Kistler  8620  accelerometer 

was  used  to  measure  the  response  at  ^  =  0°  for  various 
axial  positions  (x=0, 1.33, 2.67, 4,...  14.67, 16  ia).  The 
magnitude  of  the  displacement  response  of  the  cylinder 
at  zero  angle  for  x=0, 1/2,  and  1/4  can  be  seen  in  Fig.  5. 
The  plot  diows  where  the  resonant  frequencies  of  the 
cylinder  occur  and  the  mode  which  contributes  the  most 
to  each  resonance.  The  response  at  the  ends  (x  =  0  & 
16”)  is  generally  an  order  of  magnitude  less  than  the 
maximum  transverse  motion  at  a  particular  frequency. 
The  other  observation  that  can  be  made  from  Fig.  5  is 
that  the  accelerometer  located  at  the  cylinder’s  center 
does  not  detect  motion  at  the  (3,2),  (4,2),  and  (2,2) 
modes. 


Property" 

Cylinder 

PZT  Actuator 

Young’s  Modulus 
(Pa) 

64x10’ 

63x10’ 

Density  (Kg/m^) 

2700 

7600 

Poisson’s  Ratio 

0.3 

0.3 

Figure  5.  Cylinder  Displacement  at  ^  =  0  deg. 


The  displacement  magnitude  at  zero  angle  along  the 
length  of  the  cylinder  is  shown  in  Fig.  6.  Clearly  the 
low  frequency  resonances  correspond  to  the  axial  “one” 
modes  (i.e.  (3,1)  mode)  and  the  higher  frequency 
resonances  correspond  to  the  axial  “two”  modes  (i.e. 
(4,2)  mode).  The  other  important  observation  that  can 
made  from  the  plot  is  the  small  displacement 
occurring  at  the  ends  of  the  cylinder.  This  agrees  with 
the  theory  that  the  boundary  conditions  are  indeed 
simply-supported. 

Displacement  Magnitude  along  the  Cylinder  Length 


Frequency  jHz] 

Figure  6.  Spatial  &  Spectral  Cylinder  Displacement 

The  values  of  the  modal  indices  have  already  been 
placed  next  to  the  associated  resonance  in  Fig.  5  &  6  for 
clarity.  However,  the  circumferential  modal  index  for  a 


6 

American  Institute  of  Aeronautics  and  Astronautics 


given  frequency  can  not  be  determined  without 
analyzing  Ae  circumferential  displacement.  For  each 
of  the  resonant  frequencies  in  Fig.  5,  the 
circumferential  displacement  is  determined.  Figures  7- 
12  show  the  circumferential  and  axial  displacement 
magnitude  for  a  given  resonant  frequency.  Based  on 
the  displacement  patterns,  the  dominant  modal  index  is 
determined.  The  experimental  results  indicate  the 
cylinder  is  vibrating  similar  to  a  cylinder  having 
simply-supported  boundary  conditions. 


X  10"*  Displacement  Magnitude  Along  the  Cylinder  Length 
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Figure  7.  Cylinder  Displacement  at  910.5  Hz  (2,1) 


Figure  9.  Cylinder  Displacement  at  1533  Hz  (4,1) 


-jQ-«  Displacement  Magnitude  Along  the  Cylinder  Length 
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Figure  10.  Cylinder  Displacement  at  1819  Hz  (3,2) 


Figure  8.  Cylinder  Displacement  at  965.7  Hz  (3,1) 
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X  10**  Displacement  Magnitude  Along  the  Cylinder  Length 
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Figure  1 1.  Cylinder  Displacement  at  2039  Hz  (4,2) 
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iq*7  Displacement  Magnitude  Along  the  Cylinder  Length 
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Figure  12.  Cylinder  Displacement  at  2185  Hz  (2,2) 

The  analytical  structural  model  is  used  to  predict  the 
displacement  response  of  the  cylinder  and  is  compared 
to  die  experimental  results  in  Figures  13  and  14.  Table 
2  compares  the  natural  frequencies  of  the  analytical 
model  based  on  Love’s  equations  to  the  experimental 
results  and  also  to  the  frequencies  predicted  by  finite 
element  analysis  (FEA).  The  FEA  was  performed  with 
the  software  I-DEAS  using  30  axial  and  72 
circumferential  thin  shell  elements.  The  analytical  and 
FEA  natural  fi?equencies  are  in  agreement  to  within 
less  than  2.25%  for  all  modes.  The  experimental 
resonant  frequencies  and  the  analytic^  natural 
frequencies  are  in  excellent  agreement  for  the  (3,1), 

(4.1) ,  and  (4,2)  modes.  Less  accurate  coincidence  is 
found  for  the  (2,1),  (3,2),  and  (2,2)  modes.  Some  of 
these  differences  can  attributed  to  nine  flats 
(1.5x1.25”)  that  were  milled  on  the  outside  of  the 
cylinder  to  accomodate  for  other  actuators.  An  extra 
resonant  peak  appears  at  1328  Hz.  This  extra  peak  can 
be  attributed  to  the  non  ideal  boundary  conditions.  The 
only  resonant  frequency  that  has  edge  displacements  of 
similar  order  as  the  axial  maximum  displacement 
occurs  at  1328  Hz.  It  can  therefore  be  concluded  that 
this  is  a  resonant  frequency  of  the  mechanical  boundary 
(annular  shim).  The  experimental  and  analytic^ 
responses  of  the  cylinder  do  exhibit  similar  trends  in  the 
location  of  the  resonant  frequencies  and  in  the  influence 
of  the  actuator  location  on  the  response.  Clearly  the 
boundary  conditions  created  in  the  experiment  are  not 
perfect;  however  they  do  provide  a  good  comparison 
for  the  analytical  model.  More  importantly  for  the 
acoustic  problem,  it  is  necessary  for  the  structural 
analysis  to  be  able  to  predict  the  amplitude  of  the 
velocity  response.  By  inspecting  Figures  13  and  14  it  is 
clear  that  the  analytical  and  experimental  displacement 
magnitudes  are  very  similar  for  the  (3,1),  (2,1),  (4,1), 

(3.2) ,  and  (4,2)  modes.  This  is  of  fWdamental 
importance  in  order  to  determine  the  authority  of  PZJ 
actuators  for  the  rocket  fairing  problem. 


Mode 

Analytical 

FEA 

Exp. 

%  Error 

(3,1) 

981  (Hz) 

973 

966 

1.5 

(2,1) 

1084 

1082 

911 

19.0 

(4,1) 

1530 

1525 

1533 

0.2 

(3,2) 

2022 

2014 

1819 

11.1 

(4,2) 

2082 

2060 

2039 

2.1 

(1,1) 

2190 

2191 

- 

- 

(5,1) 

2378 

2356 

- 

- 

(5,2) 

2740 

2679 

- 

- 

(2,2) 

2798 

2792 

2185 

27.6 

(4,3) 

3013 

2946 

- 

Table  2.  Comparison  of  Natural  Frequencies 


Frequency  [Hz] 


Figure  13.  Analytical  &  Experimental  Response  at 
x  =  l/4,^  =  0” 

For  the  analytical  natural  frequencies  which  correspond 
to  (4,1)  and  (3,2)  modes.  Fig.  15  and  16  show  the 
analytical  spatial  displacement  response.  The 


Frequency  [Hz] 

Figure  14.  Analytical  &  Experimental  Response  at 
x  =  l/2,0  =  0° 


8 

American  Instimte  of  Aeronautics  and  Astronautics 


Phase  Adjusted  Cylinder  Response 


Anaular  Position.  Tbsta  idea.  I 


Figure  15.  Analytical  Response  at  1533  Hz  (4,1) 
Phase  Adjusted  Cylinder  Response 


300  0  '*'*'®*  Positic 

Angular  Position,  Theta  |deg.) 

Figure  16.  Analytical  Response  at  2021  Hz  (3,2) 

phase  of  the  displacement  response  has  been 
incorporated  into  the  plots  to  show  actual  normal 
displacements.  As  is  expected  for  each  of  these 
frequencies,  the  displacement  response  is  similar  to  the 
dominant  mode.  The  results  indicate  that  the 
impedance  model  for  PZT’s  actuating  a  SS  cylinder  can 
pr^ct  the  overall  resonant  and  modal  behavior  but, 
more  importantly  for  acoustical  applications,  is  able  to 
determine  the  magnitude  of  the  vibration. 

Future  Work 

The  structural  model  needs  to  be  further  verified  for  a 
broader  frequency  range  with  more  experimental 
information  about  the  absolute  spatial  displacement  of 
the  cylinder.  This  can  be  accomplished  by  using 
scanning  laser  vibrometiy.  The  model  also  needs  to  be 
verified  for  out  of  phase  actuation  and  more  attention 
needs  to  be  given  to  the  actual  damping  values  on  a 
mode  for  mode  basis.  The  results  from  the  impedance 


model  will  ultimately  be  used  to  calculate  the  acoustic 
response  within  the  cylinder.  Once  the  acoustic  model 
is  verified,  the  combined  structuial/analytical  model 
will  be  used  to  determine  the  authority  of  PZT  actuators 
on  a  cylinder  with  properties  similar  to  a  rocket  payload 
fairing. 

Conclusions 

An  experiment  of  a  cylinder  with  simply-supported 
boundmy  conditions  has  been  conducted  and  was 
shown  to  exhibit  dynamic  behavior  very  similar  to  what 
the  theory  predicts.  An  analytical  impedance  based 
model  of  a  PZT  actuating  on  a  SS  cylinder  was  also 
described  and  compared  with  the  experimental  results. 
The  results  indicate  that  the  imped^ce  based  model 
can  predict  the  overall  resonant  and  modal  behavior 
but,  more  importantly  for  acoustical  applications,  is 
able  to  detemtine  the  magnitude  of  the  vibratioa  An 
acoustical  formulation  of  the  pressure  field  within  the 
cylinder  has  also  been  described  based  on  the  Kirchoff- 
Helmholtz  integral.  The  work  presented  provides  the 
iititial  steps  required  in  determiiung  Ae  acoustic 
authority  and  feasibility  of  using  piezoelectric  actuators 
on  rocket  fairings. 
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